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Abstract 

In lattice Hamiltonian systems with a quartic coupling 7, a critical value 7* may exist such 
that, when 7 = 7*, the leading irrelevant operator decouples from the Hamiltonian and the leading 
nonscaling contribution to renormalization-group invariant physical quantities (evaluated in the 
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the lattice Hamiltonian of vector spin models with O(N) symmetry. As a byproduct, systematic 
asymptotic expansions for the relevant lattice massive one-loop integrals are obtained. 
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I. INTRODUCTION 



The quest for better analytical and numerical methods in the theoretical evaluation of 
measurable physical quantities, like critical exponents and amplitude ratios, is one of the 
lasting tasks of statistical field theory 

In recent years substantial progress in this field has been made by the introduction of a 
method based on the strong coupling lattice expansion of improved Hamiltonians [|T| (for a 
review, cf. Ref. 0). The essential feature of this method is the possibility of removing all 
leading nonscaling contributions to physical quantities, in the neighborhood of criticality, 
by a specific choice of a parameter in the lattice Hamiltonian (critical coupling). The 
convergence of analytical and/or numerical evaluations is therefore impressively faster than 
in any other variant form of the models under inspection belonging to the same universality 
class. 

The main limitation of this method lays in the absence of an efficient analytical technique 
for the determination of the critical parameter. Conceptual reasons for this limitation may 
be found in the impossibility of an e expansion for the value of the critical parameter, both in 
the 4 — e and in the 2 + e expansion schemes. In practice one must resort to an extrapolation 
from numerical Monte Carlo finite-size evaluations of some physical quantity, typically the 
Binder cumulant |], 0. 

Another well-known analytical approach to the study of critical lattice models is the 1/N 
expansion, which applies in particular to the physically very important class of O(N) spin 
models in three dimensions. However it is known that, in exactly three dimensions and 
for nearest-neighbor interactions, the critical parameter in the large- N limit is a negative 
number |§ . Since the critical parameter controls the large- field behavior of the interaction 
potential, a negative value would naively imply an unbounded Hamiltonian. In practice 
this would prevent a Monte Carlo simulation of the system, thus seriously jeopardizing the 
conceptual meaning of the whole approach. 

Nevertheless we still seriously believe the 1/N expansion to be, from a theoretical point 
of view, probably the most relevant expansion scheme that can be applied to any quantum 
and statistical field theory, in that there is no known obstruction to summability of the series 
expansion in powers of 1/N for the values of physical quantities. 

Therefore we decided to explore the conceptual and numerical consequences of performing 
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a systematic 1/N expansion of the critical parameter, for the class of three-dimensional O(N) 
spin models, in order to check the actual relevance of the drawbacks that we mentioned 
above. 

As a consequence of our analysis we found that the critical parameter can be formally 
computed within the expansion with no limitation related to the sign of its large- A r value 
and for space dimensionalities in the interval 2 < d < 4. In particular we found that the 
sign of the first 1/N correction is positive, and one may then expect to find a value N c such 
that the parameter itself vanishes. It is however still unclear that the predictions of the 1/N 
expansion may be extended to the region N < N c . 

As a byproduct of our analysis we obtained new more efficient expressions for the asymp- 
totic expansions of many important functions entering our calculations. We presented these 
results with some details because they might be relevant to other computations of critical 
and subcritical quantities. 

In Sec. m we introduce the lattice O(N) models and their 1/N expansion. In Sees. [II, 
IV| , and [V] we discuss on general grounds the relevant asymptotic epansions (gap equation, 



effective propagator and renomalized coupling). In Sees. |VT| and |V1J| we specialize these 
expansions to the case of the standard nearest-neighbor interaction, with the help some 
useful integral representations. A few numerical results are presented in Sec. Vlll . In Sec. 
FX] we compute the 1/N correction to the improvement parameter and finally in Sec. [X| we 
discuss the meaning and relevance of our results. 



II. THE EFFECTIVE HAMILTONIAN AND THE GRAPH EXPANSION 



Our starting point will be the usual iV-component 4 lattice Hamiltonian in d dimensions: 



(1) 



where V^0(x) is some (local) form of the lattice gradient; in the standard nearest-neighbor 
formulation V M 0(x) = <fi(x+fi) — <f>(x). 

Following Ref. H we define the rescaled couplings 



6^o 



7 



9oN' g N 

and introduce an auxiliary field a in order to eliminate the quartic term in the Hamiltonian. 
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After a trivial Gaussian integration the resulting effective Hamiltonian is 

H ee = y[TrLn/3(-V M V M + ia) - if3a + ^a 2 }. (2) 

In the limit 7 — > this Hamiltonian reduces to the usual effective large- N expression for 
the nonlinear a model. In the nearest-neighbor formulation, it is also known as the O(N) 
Heisenberg model; its large- N limit was investigated in Ref. ||. 

The saddle-point condition on the effective Hamiltonian leads to the so-called gap equa- 
tion 

9 r d d p 1 
^ + 27^0= / /n L^ ; -s , (3) 



w (27r) d p 2 + m, 



where p 2 is the Fourier transform of the lattice Laplacian operator — V^V^, which in the 
nearest-neighbor case takes the form p 2 = 2 (1 — cosp^). 

The gap equation allows for the elimination of ft in favor of the new parameter mo (large- 
iV inverse correlation length) in the Feynman graph expansion. 

In the large- iV limit criticality corresponds to the vanishing of Wq, and the criticality 
condition may then take the form 

r d d p 1 

corresponding to a finite value of (3 C for all d > 2. 

The graph expansion for this model, in the formulation based on the effective Hamil- 
tonian, requires defining the (bare) propagator A for the effective field a; by standard 
manipulations one obtains 

D(k,m , 1 ) = A- 1 (k,m n) = l:[ f \ d 2 j 2 1 == j= ; j + 7- (4) 

The above approach is quite general, and it leads to a systematic 1/N expansion of the 
correlation functions and of the physical quantities for arbitrary values of thq and 7. 

However we want to focus our attention on the critical domain, and in particular we want 
to evaluate the coupling 7* such that the first nontrivial corrections to scaling turn out to 
vanish in the computation of physical quantities in the scaling region. 

To this purpose it is convenient to parametrize the cutoff dependence of correlation func- 
tions and renormalized couplings in terms of the lattice spacing a, which can be made to 
appear explicitly in calculations by a rescaling of the coupling and momentum dependence. 
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Around criticality the dependence on a is not analytic, and as a consequence we need 
asymptotic expansions in order to identify the scaling and leading nonscaling contributions 
to any computable quantity. 

The basic technique for asymptotic expansions in powers of m^a is described in Ref. ||; 
here we shall discuss its applications to the cases of interest for the present paper. 

We only recall that, in order to regularize the generic lattice integral 

r d d p 

I(k;m a) = J j^-^F(t,m a,p), 
where k is any collection of external momenta, we can make use of the formal identity 



where 



I(k; m a) = I\ at {k; m a) + I con (k; m a), 



/ n d d p f + °° d d p 

m a, p) - J a d j^-yT m F (h, m a, pa) 



and 



/ dp 
m a,pa) — T^ v ^F(k; m a, pa)]. 
-oo (A) 

The T = T^ IR ) + T( uv ) operation amounts to a Taylor series expansion of the integrand 
in powers of m^a, and T^ uv ^ corresponding respectively to the IR and UV singular 

terms in the expansion. 

It is possible to prove that the expansions of 7i at (the "lattice contribution") and I con 
(the "continuum contribution") are individually and fully regular; the nonanaliticity of the 
expansion is factored out in the a d term multiplying the continuum contribution. 

III. ASYMPTOTIC EXPANSION OF THE GAP EQUATION 

For the purposes of the present paper and in order to show an explicit example of the 
asymptotic expansion procedure let us consider the gap equation in the large- iV limit. 

Let us assume a generic lattice Laplacian such that p 2 ~ p 2 + cp 4 + 0(p 6 ), where p 2n = 
YlnPfT' i n the nearest-neighbor version c = — ^. 

From the previously derived results we obtain the following relationship: 

" w d d p 1 



2„2 



(3 C - (3 = m a 



n (2-K) d p 2 (p 2 + mla 2 ) 



+ 2 7 



Let us now define 



G(p, m a) = 



m^a 2 



p 2 (p 2 + m^a 2 ) 

and perform the relevant expansions up to the first few nontrivial terms: 



,2^2 



4^4 



T G(p, mod) 



m^a rrina 



(p 2 ) 2 (p 2 ) 3 



G(pa,m a) ps — 



T G(pa, mod) 



1 1 
1 



p 2 p 2 + UIq 



cp 



rric 



(p 2 ) 2 (p 2 ) 



2 4 

— m cp 



(p 2 ) 2 (p 2 + 771q) 2 

2 



3m 2 



(p2)3 (p2)4 



+ 0(a 2 ), 
+ O(m 6 ,a 2 ). 



By grouping together the IR singular terms we therefore obtain the lattice contributions: 



2m 2 a 2 ( 1 - 7o ) - m 4 a 4 <5 + O(m b a b ), 



where we defined the following numerical constants: 



1 



7o = 



So 



n d d p 1 



+oo id 



dp 



w d d p 1 



> {^Yip 2 ) 2 

d d p ( 1 



3c 



p 



(5) 



.A^YifY J-oo (2nY\(p 2 Y ~~(p 2 ' 4 

On the other side, by grouping together the UV singular terms, we obtain the following 
continuum contributions: 

" + °° d d p 



(2tt) c 



d-2 



1 



p 2 + m\ p 2 



+ a d cp 4 



1 



2 

mi, 
2- 



{p 2 + Wq) 2 (p 2 ) 2 (p 2 ) 



+ 0(a d+2 ) . 



We can perform the continuum integrals by standard dimensional regularization tech- 
niques. We then finally find 

C - (3 « -6q (m a) d ~ 2 + 2( 7 - lo )rr%a 2 " |c&c>mj{a d - 5 m 4 a 4 + O((m a) d+2 ), (6) 



where 



6 - r(1 "^ 
o = 



(47r)i 



From the asymptotic expansion of the gap equation we immediately learn the following 
lesson: it is possible to choose for the quartic coupling 7 the special value 70 such that the 
first nonscaling contribution to the large- N saddle-point condition vanishes. 

By this procedure we have identified the large- N critical coupling for these versions of the 
model. In fact any change in the form of the local interaction, reflecting itself in the detailed 
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FIG. 1: 70 vs. c. 

p dependence of the bare lattice massless propagator without changing its singular part, leads 
to a finite change in the numerical value of 70 without affecting its formal representation. 

It is especially interesting to consider the class of models characterized by next-to-nearest 
neighbor interactions, and whose propagator of the fundamental excitations is obtained in 
the form 

d 

P 2 = Y1 ( 6c + 2) - 8 ( c + D cos ^ + 2 ( c + T2) cos 2 Pn =P 2 + {c+ ±)p\ 
,t=i 

where 

d 

Note that c = — i corresponds to the standard nearest-neighbor interaction, while c = 
corresponds to the 0(a 2 ) Symanzik tree- improved version of O(N) models M. 

In three dimensions we have numerically explored the range — ^ < c < 0: our results are 
presented in Figs. |1] and |2|. 

Let us notice in particular that the choice c = —0.03332110... corresponds to vanishing 
70, and it is therefore, at least in the large-iV limit, an alternative version of a spin model 
where the leading corrections to scaling are automatically made to vanish. In turn when 
c = we obtain 70 = 0.003328210..., which implies a small but nonvanishing leading scaling 
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FIG. 2: 5 vs. c. 

violation, and do = 0.002181406... . Compared to the standard this version has 

however the advantage of being numerically testable also by Monte Carlo methods, since 
7o > 0. 

IV. ASYMPTOTIC EXPANSION OF THE EFFECTIVE PROPAGATOR 

In order to perform an asymptotic expansion of the 0(1/N) contributions to physical 
quantities we must compute one-loop graphs involving the effective propagator A(k, m ,7). 

We therefore need to evaluate the asymptotic expansion of A or, more conveniently, of 
D(k, mo, 7). It is easy to recognize, from the definition of D and from the general rule of 
the asymptotic expansion, that in general we may express the result in the form 

00 

D(k,m a,j) = Y} A n{k){moa) 2n + B n (k)(m a) 2n+d - 2 ], (7) 

n=0 

where A n (k) have the form of lattice contributions (and only A Q depends on 7), while 
B n {k) are continuum contributions that can be analytically computed, e.g., in dimensional 
regularization. 

We shall not repeat the derivation (some details can be found in Ref. ||), but only quote 
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the relevant results: 
A x {k) 



+ 7, 



i r d d P 1 i 

r ^ i i /"+ 00 d d P i i 

" y_„ (2n) d (p 2 ) 2 (p + fc)2 + 7^ (27r)V) 2 F' 

J /"+oo 
-T) 2 ) 2 " loo 



" d d p 1 



1 1 



d d p 



1 



) rf (p 2 ) 3 (p + *0 2 2(p2)2 ((p+A;)2) 



3c 



p 



(p2)4 fc 2 (p2)2 



+ — — C (A;) + - 



1 



+ 



1 



( 27 r)^(p2)3^ 
1 



where we defined 



B 1 (k) = -b 



2 V(p 2 ) 2 ((p + k) 2 ) 2 J (fc 2 ) 2 ' 
! C (k) 



TO, 



2 ' 



2 A; 2 



(8) 



Trivial manipulations allow to express A n (k) as pure lattice integrals: 



A 1 (k) 
A 2 (k) 



w d d p 



1 



k 2 (p + kf 



+ 



/7T 
-7T 



ff d d p 

n d d p 
(2ir) d 



+ 2 



1 



To 
k 2 ' 



1 



{ff (p~+k)2 (p2)3 £2 (p2) 
11 1 111 



2 C (fc) 



+ ^-27oC (A;) 



(9) 



1 



- 2- 



7o 

\k 2 y 



_2 {p 2 ) 2 ((p + k) 2 ) 2 2{p 2 ) 2 (k 2 ) 2 2 ((p + k) 2 ) 2 (k 2 ) 2 _ 
From our general considerations it should be by now clear that we shall also need a 
different expansion of D, homogeneous in powers of to and k. 
Without delving into the details, we find that 



1 

D(ka, TO a, 7) ~ -a 



d-4 



+00 Jd 



dp 



1 



1 



+ (7-7o) + 0(a d - 2 ). 



(2n) d p 2 + ml(p + k) 2 + m 2 , 

Notice that the integral can be analytically computed in all dimensions, and the final 
result is 

. > r.' r.' <2 



D(ka, TO a, 7) ~ d ( -^k 2 a 2 + m^a 2 



' ( : -till) i - ---n », 



2 1 2 1 2 7 e 



(10) 



where do = |(1 — f)&o, 2 -Pi is the Gauss hypergeometric function and £ = ^/l + 4toq//c 2 . 

Let us in general denote by the label "c" the quantities occurring in the leading order in 
the homogeneous (continuum) expansion of D, and in particular 7 C = 7 — j . 
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Following Ref. ||, we can exploit identities between hypergeometric functions to recast 
the above result into the form 

-d-3/r.„\d-4 . , fro fd — 1 d 1 \ d _ 2 



D c = a H d -\ka) d -* + 7c + *F X ( — , 1, - , 1 - ~ 2 ) (m a) d " 2 , ( I I ) 

where 

ir(f -i) 2 r(2-f) 



2 (47r)lr(d-2) 
For d = 3, Eq. (0) reduces to 

1 1 jfc 

-Dc = T7 — r + 7c — 75 — r arcc °t 



16 ak 8nak 2m 

Eq. ( PD is especially appropriate for the asymptotic expansion of -D c , because all the 
nonanalytic dependence on mod is explicitly factored out in the last term. In particular 
(after a rescaling ka — > fc) we obtain the following behaviors: 

A 0c = «o(A: 2 )^ 2 + 7c 4* = -^ r 4 T!^\ ao(fc 2 )^ 2 -". 

n! — n) 

Finally notice also that the zero-momentum value of D is related to the derivative of the 
gap equation with respect to the mass, and we can obtain the relationship 

D = D(0, m a, 7) « d (m a) a! " 4 + 7c - ~^cb (m a) d ~ 2 - 5 m 2 a z \ (12) 

o 

and as a consequence 

L>oc = d (m a) d ~ 4 + 7 C . 

V. ASYMPTOTIC EXPANSION OF THE RENORMALIZED COUPLING 

We now recall from the literature the expression of the 0(1/N) contribution to the 
(unrenormalized) self-energy of the fundamental quanta: 

Ei(p,m ) = E la (p,m ) + E lb (p,m ), 

where 

f d d k A(t,m , 7 ) 

1 d d A; 9 

Sife(m ) = -A(0,m ,7) J 7^p A ( fc ' m °' T) ^I^^" 1 ^' m °' 7 )' 
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Again we might perform an asymptotic expansion of this expression, on the lines traced 
in Ref. |J, recovering in the scaling limit the (unrenormalized) continuum contribution, and 
in principle we might evaluate the first nonleading contribution. 

However for our purposes it is much more convenient to work directly with quantities cho- 
sen in such a way that all renormalization effects are automatically removed, i.e., quantities 
whose scaling limit is a finite, renormalization- group invariant, amplitude. 

The simplest such object is the so-called "renormalized coupling" g r , whose continuum 
(scaling) value g* has been computed to 0(1/ N) in Ref. [§]. 

Actually the formal expression derived in Ref. @ is correct also for the lattice ver- 
sions of the model, when continuum propagators are replaced by their lattice counterparts. 
Parametrizing the result in terms of the renormalized mass m we therefore obtain 

1 , I 



g r (m,i) = m A(0,m,7) 



l + -^ 1) (m, 7 )+0 



N 2 



(13) 



where gi is given by 



£ ; (m, 7) = A(0.m, 7) ( S la (0, m) + S 16 (m) - m 



+ 



dm 2 ' 

2 9S la (0,m) d£ 16 (m) 



and we defined 



dm 2 



T(m ) 



dm 2 
d d k 



dp 2 



dp 2 

- 2A~ 1 (0,m,7)T(m 



(14) 



A(fc,m ,7) 



_ w (2vr) d V A; 2 + mg 

It is now a matter of trivial algebraic manipulations to show that the expression for gr 
can be cast into the form 
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(D = JL 



d d fc [1 (1 d 



D J (27r) d [D \2dm 2 [dm 



3D 



+ 



2D, 







1 



2D 2 



2D, 



k 2 + m 2 
n 2 



Dr 



(k 2 + m 



2^2 



2D + m- 



dP 
dm 2 



C 







(15) 



9m 2 fc 2 + m 2 

In the asymptotic expansion of we may again identify a continuum and a lattice 
contribution. 

Renormalization-group theory insures us about the expected properties of these contribu- 
tions. In particular the leading continuum term should be finite and should therefore require 
no UV counterterms. Its evaluation will however allow us to identify the IR counterterms 
needed in order to regularize the lattice term. 
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The lattice contribution in turn should not affect the leading (scaling) order, since g* is 
an invariant amplitude which should not depend on the detailed form of the Hamiltonian. 
In turn, trivial power counting in a shows us that the first correction to scaling is generated 
by the leading lattice contribution, which is 0((ma) A ~ d ). 

We can prove the following identities, corresponding to similar results of Ref . [|Q : 
dD r 2 



dm 2 k 2 + 4m 2 
dD 0c 



[(d _ 3)D C - D 0c + (4 - d) 7c 



dm 2 

C c 



i_ 2 ) D ° c ~ 7 ' 
2 / ™2 



4 k 2 



- 1 



[k 2 + m 2 ) 2 \dk 2 + m 2 
As a consequence we obtain in leading order the following continuum contribution, de- 
pending on the single (dimensionless) variable x = / y c (ma) 4 ~ d : 



9r,con( X ) 



Oc 



+oo 



d d k 
J2kY 



d-3 



+ 



7c 



4-d 



+ 



Oc 



3m 



k 2 + 4m 2 D c k 2 + Am 2 D c (k 2 + m 2 )(k 2 + Am 2 ) 



+ 2 



(d-3)(d-5) 



+ 



+ 



7c(4 - d) 

2 7c (4-d)m 2 
lT c d{k 2 + rri 



2(d-5) 



+ 



(F + 4m 2 ) 2 2(F + 4m 2 )(fc 2 + m 2 



+ 



1 - 



(P + 4m 2 ) 2 

^ 1 

d/ (fc 2 + m 2 ) 2 



2A3 



Oc 



m 



D c k 2 + m 2 



6(d- 5) 



d- 8 



(16) 



_(P + 4m 2 ) 2 2(fc 2 + 4m 2 )(F + m 2 ) (P + m 2 ) 2 _ 
which can be shown to correspond exactly to the result presented in Ref. HQ. 

In particular the fixed-point value is obtained by setting x = 0, corresponding to 7 C = 0, 
that is the condition for the removal of the leading nonscaling behavior in the large- N limit. 

From the above result we may immediately read off the structure of counterterms, just 
by taking the power series expansion in powers of m 2 , and in particular by replacing D c 
with A Qc and D Qc with do(ma) d_4 . The resulting (singular) expression is 

2 



[ma 



+oo 
A—d r Ad 



do 



d d k 
(2vr) d 



-2 



d-3 j^A_ 

k 2 



d 



A 0c k 2 



(d-3)(d-5) 7 . 
(k 2 ) 2 A 



Oc 



d-4 



1\ 4-d 



d; (A; 2 )^ 



Let us now compute the lattice contribution. By applying the already defined asymptotic 
expansions we obtain: 



dD 



+ 



2Dn 



dm 2 k 2 + m 2 



A 1 + 2^-bo{*C + 



(k 



2\2 



m 



d-2 



+ 2 A 2 - 



7c 



2^2 



5o 
k 2 



m 2 + 0(m d ) 
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and 



1 d 



2 dm 2 



dD 



+ 



2Dn 



dm 2 k 2 + m 2 



Dn 



A 2 -|-2 7c (C + 



(k 2 y 



(k 2 + m 2 ) 2 
+ 0(m d - 2 ) 



2D + m 



2 dDo 
dm 2 



C 



Notice that the most singular contributions, 0(m d ~ A ), have been removed in the above 
combinations. 

Let us now find out the leading IR singularities of the terms appearing in the integral 
representing g^iir We mus t onr y notice that, for all A n (k), the singular behavior when k — > 
is determined by the corresponding behavior of A nc (k), with corrections whose singularity 
is depressed by a factor k 2 . As a consequence we obtain in the IR limit 



-^U 1 + 2^) --^ 



A 



k 2 



k 2 



d-3 + (4-d)-^- 



(17) 



1 

A 



(k 2 y 



(d - 3)(d - 5) + d - 4 - - (4 - 



7c 



A 



Oc 



These singularities are perfectly matched by the terms coming from the expansion 
of the continuum contribution, as expected. As a consequence we are able to write down an 
exact, finite representation of the leading lattice contribution to taking the form 

\4-d 



&( m >7c) = 



ma 



do 



(18) 



where 



d k_ 

Y 

oo 

d d k 



1 



A 2 - 5 S-2 lc \ C -} 



(k 2 y 



{27i) d (k 2 ) 2 



(d - 3)(d - 5) + d - 4 - - (4 - d) 



lc 



A 



Oc 



d - 3 + (4 - d) 



lc 



1 2 



A 



Oc 



(19) 



It is worth observing that ^con(^) showed a nonanalyticity whose leading dependence was 
proportional to xhix. Correspondingly Sg^^c) has a nonanalytic 7 c ln7 c dependence. The 
coefficients of these singularities match properly in order to reproduce an overall dependence 
proportional to r y c (ma) 4 ~ d In ma, as expected from general renormalization group arguments 
because of the anomalous dimension of the leading irrelevant operator. 
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VI. INTEGRAL REPRESENTATION FOR THE GAP EQUATION 



For the purpose of actual numerical calculations one must find an efficient way of perform- 
ing lattice momentum integrals. In practice this may be obtained by resorting to parametric 
(Feynman and Schwinger) representations of the lattice propagators. These representations 
are especially useful in the case corresponding to the standard nearest-neighbor Hamiltonian. 

Let us first consider the integral appearing in the gap equation 



n d d p 1 



Introducing Schwinger's proper time representation we obtain 

roc ptt id^ poo 

Z(m 2 )= dae~ a < ^_^ e -2-E M (i-cos Pfl) = / da e ~*< \ e -2<* Io ( 2a )] , (20) 
Jo J-* ( 27r ) Jo 

where Iq is the standard modified Bessel function, admitting for large values of its argument 
the following asymptotic expansion: 

/nU (4™)*^ nir(i-n) (4a)»- 



V (2a) w— —rJ2 



By a proper change of variables the numerical evaluation of the above integral is now possible 
even in the small regime. 

When d — 3 it is also possible to obtain an analytical expression for 2(0) ||, which was 
first derived in Ref. |J: 

& = H(0) = ^±V(x), (21) 

where K is the complete elliptic integral of the first kind and x = (2 — v / 3)(v / 3 — \/2), and 
the resulting numerical value is (3 C = 0.252731009858663... 

We can compute H(mg) as an asymptotic expansion around = 0: 

oo 

5K) = a n< n + hml™. (22) 

n=0 

It is conceivable that when d = 3 the coefficients a n of the expansion can be computed 
analytically in terms of elliptic integrals, but we contented ourselves with a numerical cal- 
culation of the "lattice terms" appearing in the analytic part of the expansion, while it is 
always possible to obtain closed form expressions for the coefficients b n of the nonanalytic 
part. 
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We can obtain explicit expressions for the coefficients in the expansion fl22"|) by subtracting 
a proper number of terms of the asymptotic expansion of the Bessel function raised to the 
power d. Let us label the coefficients of this expansion according to the equation 



-2a 



/o(2a) 



a 



(4vra)^ n=0 

the coefficients can be computed recursively from the equations 

co{d) = 1, 

^(d) = - ^(*d - n + fc) c n _ k . 

k=l v 2 / 

Let us add to the integrand in Eq. (f20D the formally vanishing term 
1 



(23) 



(24) 



OO / D\<n 71—1 

— am, v 



£ 



-^c m (d)«- m + -l- r ^c n (d)a-™ £ 
(4vra)2 n=Q 



—a m, 



2\m 

oJ 



m! 



(47Ta)- n=Q ■■• m=0 uy n=0 m= n +l 

where we have only interchanged the order of the summations in the two contributions. It is 
now possible to group the first contribution with the original integrand and recognize that 
the resulting combination defines an analytic function of tUq, since each coefficient 

n-l 



a,, 



a) 



da — 
o n\ 



-2a 



/o(2a) 



1 



-j^c m (d) 



a 



(25) 



(4vra)^ m=0 

in the corresponding power series expansion is a finite, UV and IR regulated integral for 
any d in the range 2 < d < 4. In turn, the integration of the second contribution may be 
represented, after trivial resummations and rescalings, in the form 



(4vr 



1 OO 



n=0 



where 



k n (d) 



dxx n 2 



m=0 



are finite UV and IR regulated integrals in the range 2 < d < 4; integrating by parts 
we obtain a recursive equation for k n , which can be solved immediately, obtaining k n = 



T(l — | —n), and consequently 



7T-7T c "( rf ) r ( 1 - o _n ) ' 
(47r) 2 V z 



(26) 
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It is now trivial to set d = 3 in Eq. (|24|) to obtain an explicit recursive expression for b n . 

To compute numerical values of a n in d = 3, it is not practical to use Eq. (pop, given 
the slow convergence of the integration for large a. We found it more convenient to split 
the a integration in Eq. ( [ZUD at a = 1. For a < 1 we can expand in powers of m§ under 
the integration sign and integrate term by term. For a > 1, we subtract n + / terms of the 
expansion (f2"3"D, expand the integrand in powers of j«q up to 0(rriQ n ), and integrate term 
by term; while the resulting integrals converge for any / > 0, it is useful to set / > 3 to 
ensure fast convergence of the integration for large a. The integral of the regulator can be 
computed analytically in terms of the incomplete Gamma function, and the singular part of 



its asymptotic expansion reproduces the singular part of the expansion (22). 

Collecting all numerical and analytical results, we obtain for d = 3 the following asymp- 
totic expansion: 

E(m 2 ) = p e - ^-niQ - 0.012164158583022 m\ + -j- m 3 + 0.00083776240606293 

47T 32 7T 

1 1 281 

-m\ - 0.000066743211781194 m 6 + m\ + 5.7884488124445 x 10~ 6 m^ 



2560 7T u u 430080 7T 

' ! -ml - 5.320777180475 x 10" 7 m^° + 1^ ml 1 



655360 7T u u 403701760 7T 

70289 

+ 5.107544566044 x 10- 8 mJ 2 - 19377684480 ~ m « + °K 4 )- (27) 

which gives an error smaller than 10~ 14 in the range < m\ < 0.1. 

It is immediate to extract from this expansion the 3-d large-iV nearest-neighbor model 
values 

7o w -0.0060820792915113, 5 w 0.00083776240606293. 

The asymptotic expansion of H^(rrio), the nth. derivative of S(mg) with respect to m^, 
is easily obtained from the above expression. There is an obvious precision loss, but in the 
above-mentioned range the error in the second derivative is still smaller than 10 -10 . 

In the following, we will need to compute E:*™) in a fast and accurate way, for generic 
values of m^; to this purpose, we tabulated H for values of m\ on a uniform grid with step 
h = 10~ 3 and compute ^ n \ml) by n+4-point Lagrange interpolation. 
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VII. INTEGRAL REPRESENTATIONS FOR THE EFFECTIVE PROPAGATOR 



Let us now recall from Ref. the following basic result: 

D(k, m , 0) = i/ d* £ g i K + 2£ , 1 1 _ z>mtf 



where = y 1 — se(1 — 

Using Schwinger representation we then obtain: 

1 /-l POO PIT Jd n 

D(k } m ,0) = - dx daae- am ° J!_^ e -2«E M (i-*M «*,,.) 

2 JO JO J — 7T ( 27r ) 

= i / cfe / c/aae- am o]T[e- 2 %(2a2 M )j. (28) 
2 Jo Jo „ L 

Alternatively we might use the Schwinger representation directly and obtain 

-I /*oo /*oo r 

D(k,m ,0) = - ds dt e - {s+t) < TT e~ 2{s+t) h{2^s 2 + t 2 + 2st cos k„) , 
2 Jo Jo „ L J 

that can be reduced to the previous one by the variable change s = x a, t = (1 — x) a. 

The direct numerical evaluation of Eq. (p8|) , and especially of its derivatives with respect 
to Mq, in d = 3 is difficult, particularly for small values of k or tuq; the convergence can 
be improved dramatically by adding and subtracting a symmetric combination of Bessel 
functions: 

D(k,m o ,0) = - [ dx f da«e- £,(6+m » ) [TT/ (2«z (I ) - ll(2az) 
2 Jo Jo L 

where z = + Z2 + 23); the subtracted integral and its first few derivatives with respect 
to wIq can now be computed accurately by Gauss-Legendre integration on a small grid; the 
integration in x of H' and of its derivatives is also easy, once a few singular terms of Eq. ( |22| ) 
have been subtracted. 

In order to identify the lattice contributions to the asymptotic expansion of D we must 
expand the integrand in powers of m§. The IR singularities turn into unsuppressed positive 
powers of cu, present in the large-a regime when z M — > 1. These singularities become worse 
and worse with increasing powers of m§. In order to classify them according to their degree 



17 



we need to consider the homogeneous expansion of the integrand in powers of a 1 and of 

Recalling once more the asymptotic expansion of the Bessel function, we can write down 
the homogeneous expansion in the form: 

,-ax(l-x)k 2 00 



JJe- 2Q / (2a^) 



(Aira) 



C n [ d, ax(l - x)k 2 u ] a n , 



n=0 



where in turn 



2n 



C n [d, ax(l -x)k 2 ] = ^2j nm (d)[ax(l - x)k 



2lm 



(30) 



(31) 



m=0 



and j nm (d) may show a dependence on the ratios k 2p /(k 2 ) p . 

Repeating the procedure developed in the previous section we may now obtain the fol- 
lowing decomposition: 



D(k, m , 0) = J2 Rn{k)m 2 n + D s {k, m ), 



(32) 



n=0 



where 
R n (k) 



2n\ 



1 POO 

dx J daa n+1 
(] Jo 



-ax(l-x)k 2 n 



(Ana 



a 



(33) 



H \ ~" ~V m=0 

is the part of A n (k) = R n (k) + S n (k) which shows a regular dependence on k in the k — > 
limit, and 



I /•! poo -ax(l~x)k 2 00 r 

D s (k, m )E- cfa / tfa £ e—o _ £ . 

1 J ° J ° (47ra) 2 n=0 L m=0 



, —am 



2\m 

oJ 



ml 



C n a 



l-n 



(34) 



(47ra). n=0 

is a calculable expression including all the singular dependence on k and admitting an 
asymptotic expansion in the form 



D s (k,m ) « ^[S n (A;)m^ + B n (k)m 2 



2n+d-2n 



n=0 



In order to compute D s (k,m ) let us notice that it may be also expressed in the form 

00 2n 

D s (k, mo) = ^ ^ 7nm(^)4 m (mo) ( 35 ) 



where 



4m( m o) 



2nm 



n=0 m=0 



00 g—ax(l—x)k 2 

da- 



(Ana) 1 - 



e - am o 



n-l 

E 

m=0 



(-omg) 
m! 



2\m 



ft 



1+m— n 



(36) 
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enjoys the property 



dl, 



-Ln—l.m- 



Therefore we only need to evaluate 

I 0m (m 2 ) = - / dx[x(l - x)k 2 ] m / dae""' 1 -* 2 -' ^— ; d 
Jo Jo 



(47T) 



2 



1 T(2-|+m) f 1 , [x(l-x)ib 
/ ax 



2lm 



2 



2 (47r)2 Jo [m 2 , + x(l - x)k 2 ] 2 ~2~ 

= l - r(2 ~^ ( % 2 r( 9 Y T dx (37) 

2 (4tt)I V9jfeV Jo [m 2 + x(l-x)k 2 ] 2 -? 

A trivial comparison with previous results shows that the last integral is directly related 
to the continuum propagator by the replacement k 2 — > k, and therefore, setting £ 2 = 
1 + 4ml/ k 2 we obtain 



hm{ml) = (-k 2 ) m 



d 
dk 2 



aoi d -*(k 2 )t- 2 + A- 2 F/ ^,1^,1-1 W- 



(38) 



jfc2|2'"\ 2 '"'2'" I 2 / 
where we must appreciate that the above expression is naturally decomposed into an analytic 
and a nonanalytic term, which implies that we can immediately relate all coefficients S n (k) 
to the derivatives with respect to k 2 of £ d ~ 3 (k 2 )%~ 2 . 

Straightforward manipulations lead to the general relationship 

(39) 

The above expression brings into evidence a peculiar feature exhibited by the functions 
S n (k) when d = 3. In this case the arguments of the function T(l — n + p) appearing 
in the denominator are integer nonpositive numbers whenever p ^ n, and therefore the 
corresponding contributions to the sum vanish, giving 

2 " ' U"l'i-! 



S n (k) = (-l)^(P)- £ \n_J ln q (S)- (40) 

<7=0 L V2 



9= 

When d = 3 the singular contributions to A n {k) are proportional to (k 2 )~^ for all n, and 
this is basically a consequence of the vanishing dependence on £ in the analytic part of the 
continuum propagator. 

The simplest example of this procedure is offered by the integral representation of A (k, 0): 

A (k,0) = a (k 2 f 2 - 2 + Ro(k), (41) 
19 



where 



1 f 1 f°° 
Ro(k) = - I dx I da a 



-ax(l— x)k 2 



(Ana 



(42) 



This representation is worth a few observations. Notice that the explicit term might 
be simply obtained from the expression of A 0c by the replacement k 2 — > k, and therefore 
it reproduces exactly the singular behavior of A when k — > 0. Moreover R (k) has, by 
construction, a finite limit when k — > 0. It is easy to check that Ro(0) = —70, thus verifying 
explicitly that the choice 7 = 70 leads to the 0(a A ~ d ) improvement of the lattice massless 
propagator. 

By the same technique we obtain 



Ai{k) = a 



2(d-3) + (6-d)(4-d) k A (8 - d) 



k 2 



16 



(k 2 y 



16 



(43) 



where 



1 /"l ;>oo r 

Ri(k) = — / dx da a 2 "RV 2a / (2a^) 

2 Jo Jo L „ 



-aa;(l— x)fc 2 



(47TQ;) 



1 



1 + + -x(l - x)k 2 - -ax 2 (l - x) 2 k A 

16a 4 4 



(44) 



Again we notice that the explicit term reproduces the leading singularities, and moreover 
Ri(0) = So- 

We also mention the result for A 2 (k), since it will be needed in the evaluation of 5g^: 



A 2 (k) = -oq 



O2,0 . 0,2,1 

-r 2 h + a 2 ,2 

(k 2 ) 2 k 2 



(k 2 )^ 2 + R 2 (k) : 



(45) 



where 



02,0 = 4(d — 3)(d — 5), 

(d-3)(d-12) (d-3)(d-6)(d 



02,1 = 



+ 



A- 1 



a 2,2 — 



8 (jfe2)2' 

(d-16)(d-8) (d-4)(d-6)(d-8) £ 4 



512 



+ 



256 



2^2 



(d - 4)(d - 6)(d - 8) A; 6 (d - 4)(d - 6)(d - 8)(d - 10) (k 



4\2 



64 



2^3 



512 



(k 2 ) 4 



and we obtain 



-y t-L POO " ^ — OLX[\—X)k^ 

dx daa 3 Ue- 2a I {2az^-— -j-(l + r 2 ,i(k) + r 2 , 2 (k)) 

Jo Jo (Ana) 2 



R 2 (k) = - I dx I daa i 
4 ./o Jo 



ax(l—x)k 2 



(Ana) 



(46) 
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where in turn 

r 2,l(^) = 



d 1 M 

h -x(l 

16a 4 V 



x)k 2 — -ax 2 (l — x) 2 k 4 , 



+ ^-(d + 2)x(l - x)P + 1(2(P) 2 + (8 - d)*V(l - x) 2 
512 cr 64 a 64 



_ ^ax 3 (l - x) 3 (2k 6 + k A k 2 ) + -^aV(l - x) 4 (£ 4 ) 2 . 

An obvious relationship exists between the values R n (0) and the coefficients of the analytic 
part of the asymptotic expansion of the function 5' (mo). One easily finds that R n (0) = 
|(n + l)a„+i. 

It should be by now clear that a trivial generalization of the same technique allows for 
an explicit, albeit more cumbersome, evaluation of all the functions B n (k). We could verify 
that the correct expressions for Bo(k) and B\(k) are reproduced, and we computed B 2 {k) 
for a better accuracy check of our numerical estimates. The result is too cumbersome to be 
reported here. 

In practice, A , A%, and A 2 can be computed numerically in d = 3 exploiting the sub- 
traction (^9|); since 

3k 4 - (k 2 ) 2 



-6t 



n j o(2 



-ax(l-x)fc 2 _ ,\I 



(47) 



48(4vra) 3 / 2 

the subtraction is enough to regularize A\ and A 2 . Let us consider, e.g., A\\ by straightfor- 
ward manipulations we obtain 

1 



A\{k) 



dx 



da a e 



-6a 



-1 - 



J] J (2a; 

1 



Wo*) 



167r(x(l — #)& 



2^3/2 



(48) 



it is useful to overregulate the x integration by subtracting the r.h.s. of Eq. (|47j), integrated 
over a. For A 2 we follow the same procedure; the subtractions are more complicated and 
not worth writing here. The computation of Aq is similar but of course easier. 

We verified explicitly that A n , computed in this way, and B n are consistent with Eq. ([?[) 
for n < 2. 



VIII. NUMERICAL RESULTS 



It is worthwhile to present a selection of the numerical results that we obtained in d = 
for the nearest-neighbor formulation (c = — 
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-3.50 



-3.55 

0.0 0.1 0.2 0.3 0.4 

m 

FIG. 3: (solid line) and gr}on + Q^ut (dashed line) vs. m for j c = 0. 

In Figs. [|, £|, and |] we compare gr (m, 7 C ), obtained by direct evaluation of Eq. (|15|), 
with ^r, 1 con(7c(^a) 4_d ) +5 , r 1 iit( m ' 7c), obtained from Eqs. (pUf), (|18|), and (|T9|) , for three values 
of 7 C ; of special interest is the value 7 C = —70, i.e., 7 = 0, corresponding to the nonlinear o 
model. 

Sg^^c) is plotted in Fig. || of special interest is the value 

5g {1 \0) = 0.0049699... 

IX. CONSTRUCTION OF THE IMPROVED HAMILTONIAN 

As we mentioned in the Introduction, the improvement procedure aims at a systematic 
cancellation of the next-to-leading effects in the invariant amplitudes. The renormalizat ion- 
group theory insures us that a single choice of 7 exists such that this cancellation occurs in 
all amplitudes. It is therefore sufficient to find the value 7 = 7* for which the cancellation 
occurs in the renormalized coupling. 

In the context of the 1/N expansion we may assume 7* to admit an expansion in powers 
of 1/JV: 
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FIG. 4: (solid line) and <?£con + 9 r jL (dashed line) vs. m for 7 C = —70, i.e., for the nonlinear 
a model. 




FIG. 5: gl- (solid line) and (fr.con + 9 r i at (dashed line) vs. m for 7 C = 0.1. 
We have already recognized in Sec. |TTT| that 7q =70. We may therefore define 

7c = 7* - 7o ~ t^7i*- 

Substituting this result in the expression of g r and expanding in powers of m we then 
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FIG. 6: 6gW vs. 7c . 



obtain: 

g r (ma,<y*) 



1 

do 



1 + -^SL(O) + " 7i)(ma) 4 " d + O(mV) + o(± 



We then recognize that the condition for the cancellation of the first nonleading contri- 
bution is 

" n (H- ( \ /', <\ } \ I /.4,\^ /-+ 00 d d k 4(d-3) 



7i = ^i(0) = 
where A = A (fc,7 ). 



,(27r)^Vi V 2 fcV 2 U 



(27r) d (A; 



2\2 



X. CONCLUSIONS 

In the case of the three-dimensional O(N) models with standard nearest-neighbor inter- 
actions, our analytical results led us to the prediction 

7 * -0.00608207 + 0.0049699 + O (J^J . (49) 

We might now obtain an estimate for the value iV c for which 7* vanishes. The numerical 
value of the present estimate is admittedly not very promising, but, as we mentioned in the 
introduction, the essential feature of 7^ is its positive sign, suggesting that N c may exist and 
possibly be within the range of convergence of the 1/N expansion. 
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The physical interpretation of N c amounts to the statement that, for this special value 
of N, the lattice version of the nonlinear a model (7 = 0) in three dimensions shows the 
absence of leading irrelevant operators in the expansion of the Hamiltonian into scaling 
fields. 

For all values N < N c , one would get 7* > 0, and as a consequence one may proceed 
to analyze the models in their improved version, both by numerical Monte Carlo methods 
and by perturbative expansion techniques. Numerical evidence shows that this is actually 
the case for the physically interesting cases N < 3. It might be interesting to perform a 
numerical study with the purpose of estimating N c . A naive extrapolation from the known 
numerical values of 7* (AT), 

7*(1) = 0.0159(3), 7*(2) = 0.0078(2), 7 *(3) = 0.0043(4), 7* (4) = 0.0021(7), 



obtained from Refs. M, [IU|, [IT], and |T2|] respectively using the formula 7 = /3^/(8XN), 
suggests N c ~ 5 -T- 6. 

Concerning the extension of the 1/N expansion itself to the region iV < N c , we must 
cautiously mention that some dramatic change in the analytical behavior of the function 
7*(iV) may certainly occur at iV = N c . It is not possible to perform a strong-coupling 
expansion of the models when 7 < 0, as one may immediately realize from an analysis of 
the gap equation. 

We would like to mention that the numerical evaluation of higher orders of the 1/N 
expansion is technically not beyond reach, along the lines traced by Ref. []13| and exploiting 
the more accurate results for the effective propagator obtained in the present paper. 

Finally, as a consequence of the discussion of the previous sections, it should be clear 
that, in the class of models with next-to- nearest neighbor interactions, it is always possible 
to find a choice of Hamiltonian parameters such that improvement becomes possible for 
arbitrary values of N. 
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